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Abstract 



< 

In recent experiments, the re-thermalization time of the mechanical resonator 
r^ ' is stated as the hmiting factor for quantum apphcations of optomechanical sys- 

|- tems. To explain the origin of this limitation, an analytical nonlinear investiga- 

tion supported by the recent successful experimental laser cooling parameters 
C . is carried out in this work. To this end, the effects of geometrical and the op- 

2 ' tical nonlinearities on the squeezing are studied and are in a good agreement 

^yi, with the experimental results. It appears that highly squeezed state are gener- 

ated where these nonlinearities are minimized and that high nonlinearities are 
limiting factors to reach the quantum ground state. 
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1. Introduction 

Squeezing is a beautiful quantum phenomenon with amazing potential appli- 
cations III, |_2| of which the most recent are connected with continuous variables 
quantum information [3|, \^ and ultrasensitive measurement of weak perturba- 
tions as the gravitational waves [a, la LZ| • Squeezed states are nonclassical states 

j^ I in which the variance of at least one of the canonical variables is reduced be- 

H ■ low the noise level of zero point fluctuations. To generate squeezed states, the 

- - - common technique consists to use an optical cavity filled with a nonlinear Kerr 

medium which is fed with an external pumping field (8|. With the recent ad- 
vances in cooling techniques for nano scale optomechanical systems, various 
setups have been designed for quantum ground state engineering of mechanical 
mirrors with highly squeezed states of hght [sl, llO|, [ill, llJ, [ij, |lj] . Indeed, with 



such technique it is now possible to obtain effective phonon number less than 
1 ISJ, |l6|, lul, [iSJ, [iQJ, |20| . The limiting factors to obtain much lower phonon 



'Corresponding author 
** Principal corresponding author 

Email address: snanaOuniv-ndere . cm (S.G. Nana Engo) 



Preprint submitted to Physica B April 8, 2013 



number are the re-thermalization time of the mechanical resonator Tth = -jpy- 
(where h is Planck's constant, Qm is the mechanical quality factor, ks the Boltz- 
mann constant and T the temperature of the support), which competes with 
the cooling, and the ubiquitous phase noise of the input laser which can create 
a discrepancy between experimental results and theoretical prediction 15|, Il6| . 



Nevertheless, a lot of theoretical studies on the subject has been carried out in 



the last decade and several proposals have been produced [21|, |2^, |23| . In Ref. 
|23j we applied the technique of back-action cooling to show that the cooling 
of the nanomechanical oscillator to its ground state is limited by the effects of 
both optical and mechanical nonlinearities. 

In this paper, by using the parameters of the experimental laser cooling of 



Ref. [15|, we extend the previous treatment to show through analytical study 
that there are the nonlinearities which limit the squeezing in optomechanics. 
The first one which depends on the geometry of the mechanical structure is 
known as the geo metrical nonlinearity derives from the nonlinear dynamics of 



the beams [2J, |25|, |26| . The second one is the optical nonlinearity which appears 
as a nonlinear phase shift [lO|. The geometrical nonlinearity which is always 
present and not negligible in nano resonators, is shown to be a limiting factor to 
reach the quantum ground state as suggested in Ref. [23| . In the same way, it is 
shown that high absolute values of the optical nonlinearity limit the squeezing 
of the output intensity. 

The paper is structured as follows. Section [5] will set the stage for explor- 
ing the dynamics of our system, deriving in particular the nonlinear Quantum 
Langevin Equations and the linearized equations of motion. Sections [3] and [H 
subsequently makes use of numerical simulations to discuss the squeezing of the 
mechanical and the optical output quadratures. Finally, we conclude with an 
outlook of possible future directions. 

2. Dynamics equations 



We consider an optomechanical resonator described by Fig. 1 of Ref. [15 1 



The dynamics equations of a mechanical oscillator coupled to a driven cavity is 



usually derived from a single-mode Hamiltonian [Ij, |23|, |2J, |27| as 



Fth 

Tmim + ^liX,n " P" ^7n = gM^m\a{t)\'^ + -— , (la) 

MXZPF 

a{t) ~ ie"' + y/^a'" , (lb) 



i(A + x)- - 

XZPF 2 



Xm and pm are the dimensionless position and momentum operators of the me- 
chanical oscillator related to their counterparts operators of the nanobeam as 

^ ^ V 2Mn,„ ^"' ^ ^zPFXm and p = -^r^Pm-, with [xm,P'm] = 2i. The pa- 
rameters Qm — V2wc^^^^, do, ^m and e™ are respectively the optomechanical 
coupling, the cavity length, the mechanical frequency and the amplitude of the 
input laser beam. The cavity decay rate and the mechanical damping of the 



mechanical oscillator are respectively represented by k and r,„ . The laser-cavity 
detuning is A = uii — ujc with uic the optical cavity mode frequency and ui the 
laser frequency. The terms Ftfi and a™ represent the Langevin force fluctuations 

and the input laser fluctuations. The terms (3" = ^f[''^"^ and gMXmO- repre- 
sent the mechanical and optical anharmonic terms. When the nanorcsonator is 
subjected to a large displacement amplitudes, it displays a striking nonlinearity 
/3" in its response. This comes about because the flexure causes the beam to 
lengthen, which at large amplitudes adds a signiflcant correction to the overall 
elastic response of the beam 26] . The optical anharmonic term is another kind 
of Kerr medium, which has a mechanical origin: the radiation pressure induces 
a coupling between the position of the doubly-clamped flexural resonator and 
the phase-intensity of the light beam, thus modifying the optical path known 
as the phase shift. 

One can derived from the set of equations Eqs.(IT])the following nonlinear 
Quantum Langevin Equations (QLEs) [23| 

(2a) 
(2b) 

(2c) 
(2d) 

By setting the time derivatives to zero in the set of nonlinear Eqs. (O, the 
stationary values of the position of the oscillator and the amplitude of the cavity 
field are 

r,5A/|-|2 |-|2 2K_Pm 

x„, = 2--—\a\ , \a\ =- — — — , 3 

ytm hwt ((A + gMX^Y + ^) 

The values of x obey the following third order algebraic equation, 

_3^2Axzp^_2^^^^2^^2)4pf^_^||p£i» ^0. (4) 

From Eqs. Q and (jl|), it appears that both a and Xm increase when the input 
laser power Pm increases. 

Using the experimental parameters of Ref. [15] at the detuning of A = ilm 
and for Pin = 1 mW, we obtain the following values of x which are in the range 
of those obtained experimentally in Refs. [ii,[l3]: 1-28 x 10"^^, -1.09 x 10"^ -f 



Xm = ! 


^7nPm 






Pm = -ftmXm " ^mPm + gMO^Ci + /3" + Fth 


a — 


i{A + gMXm) - 2 


a - ^£" + V^a™ 


dt = 


-i(A -1- gMXm) - - 


2. 


at -f ie* + y^a"^ 



7.43 X 10"i"i, -1.09 X 10"^ - 7.43 x 10-i°i. The first solution, which is real 
and small, corresponds to the stable regime of the mechanical resonator, while 
the two conjugate others, which have the same module {\x\ « 1.09 x 10~*), 
correspond to the unstable regime. 

For \a\ ^ 1 (satisfied in Ref. [la])i the above QLEs can be linearized 
by expanding the operators around their steady states: Xm — Xm + ^Xm and 
a = a + 5a. By introducing the vector of quadrature fluctuations u{t) = 
{5xra{t),Spm{t),5I{t),5Lp{t)y and the vcctor of noiscs n(i) = (0,Fth(t), y^5P"-{t), y^5ip'"{t)y 



where 51 = {6a' + Sa), Sip — i((5a^ — Sa) are the intracavity quadratures of the 
intensity and the phase , and the corresponding hermitian input noise operators 
SP" , Sip''"' , the hnearized dynamics of the system can be written in a compact 
form 

iiit) ^Au{t)+Ti{t), (5a) 
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The higher order of fluctuations are safely neglected. The linearized QLEs show 
that the mechanical mode is coupled to the cavity mode quadrature fluctuations 
by the effective optomechanical coupling G = gAf l^l, which can be made large 

by increasing the input laser power Pi„. /3 == ^^^f^^" and A = A + gmXm 



denote the dimensionless geometrical nonlinearity and the effective detuning. 
The range values of the geometrical and optical nonlinearities are given in the 
Table [TJ One remarks that /3 and "q increase when Xm increases and they reach 
their maximum values at the detuning A w $!,„ . As expected in the Table [TJ 
the optical and the mechanical effects are respectively high ly p ronounced at the 
optical (A « 0) and the mechanical (A « ^m) resonances [12] ■ This leads us to 
investigate the squeezing at this particular sidebands. 



Detuning A 


Mean displacement of 
the nanobcam xim) 


Range of values of nonlinearities 





2.77 X fO-" 


7;e [2.54 X 10"^; 6.79 X 10"^] 


7.42 X f 0"^" 


/3e [7.87 X 10-^5.72x lO""] 


^ '?7l 


1.27 X 10-^^ 


?7e [1.17 X 10~^1] 


1.09 X 10-« 


/3e [1.66 X 10-10; 1.22] 



Table 1 : The range of values of the optical nonlinearity r) and the geometrical nonlinearity 
at the detuning A = and A = f2„i respectively, using the parameters of Ref . [isll . 



3. Squeezing of the mechanical quadratures 

The dynamics of mechanical fluctuations is obtained by writing Eqs. ([5]) in 
the Fourier space, 

B{9)u{9)+n{n)^Q, (6a) 



where 



B{Vt) = 



f in rim \ 

rimW-i) (in^Tm) G 

(J^-f) -A 

V G A («^-f)/ 



(6b) 



Solving the matrix equation straightforwardly, we obtain the solution for the 
mechanical displacement operator to be 



Xg/f (f2)fem(f^) = aiGVlmV^. ( A^ + - - w^ + ikQ. 



-ASif'" + {-in + -] sr 



^mFth, 



where 



and 



ai 



A2 + -' - n^ 

4 



-\ -1 



'n' 



Xeff{n) = {nlff-n^-znr,ff)-\ 



(7) 



(8) 



(9) 



is the effective susceptibility of the oscillator with the effective resonance fre- 
quency and damping rate given by 



A 



nifJn)^ni^\i + a,G'— A' + - -n' -p 



K' 



l^m V 4 

2f 



(10) 



re//(r!) = r,„ - aiG'Vt„,AK. (11) 

By using the correlation functions of the noise sources for a coherent beam in 
the frequency domain, the oscillator position and the momentum variances are 
defined by. 



1 r^ 



(12) 
(13) 



where the position noise spectrum is given by 



^2r.2 '^ I a2 r>2 '^^ 






2r„if7,„fi 1 + coth 



2kBT 



(14) 



hn 



^^2^ (sat- 



At the quasi resonant frequency {fl ~ flm) and for coth \ 2k t I '^ nn 
isfied with experimental parameters used) , the exact solutions of integrals ([1 
and (fO|) are given by 



{SpD 






'^^eff^lff 



02, 



4r 



e// 



-02, 



(15) 
(16) 



where 



^2= -^} - , " . 2 ^ + 4r™ 1 + ^7^). (17) 



4nH7 i ' sF 



ftfi^ 



/ A2 
These position and momentum variances should satisfy the Heisenberg relation, 



{6xl){Spl)>\^[x„.,Pm]\\ (18) 

that is, 

{Sxl){5pl)>l. (19) 



There is no condition on the individual quadratures of relation (|19p . How- 
ever, for a standard quantum limit (SQL) of 1, the coherent states must satisfy 
(fe^„) = {Spm) = 1- When one variance is below the SQL, i.e., (fe^) < 1 or 
{Sp^) < 1, the corresponding quadrature is said to be squeezed. Generally, the 
squeezed states are characterized by the asymmetry between its quadratures 
which is mostly introduced by the nonlinear effects. According to Eqs. (I10[) 
and (jlSp . only (Sx"^) depends on the geometrical nonlinearity through the term 

nf^- Contrariwise, the value of {Sp'^) is obtained by using experimental pa- 
rameters of Ref. [l5| in Eq. pBj) and it appears to be squeezed up to about 37% 
(see FigH]). Therefore, the position variance is deduced from the mean energy 
of the nanoresonator in the steady state, 

E = ^{{Sxl,) + {SpD) ^ h^^ (^neff + ^j , (20) 

by substituting the effective phonon number with the experimental value n^ff = 
0.85 ±0.08 of Ref. [l5|. The value obtained is (Sxf^) w 4.44 which is unsqueezed 
(see Fig. HI). 

Indeed, the ratio ^yf^ increases when /3 increases (for the high mechani- 
se// 

cal displacements) and rises up the position variance {Sx'^). So, as the re- 



15| 



thermalization time of the mechanical resonator or the decoherence time [16 
the geometrical nonlinearity limits the squeezing and some quantum effects 
In fact, /3 depends on the bending moment of the resonator and takes into ac- 
count its internal vibrations. These internal vibrations increase for large bending 
and contribute to the re-thermalization of the resonator at the low temperatures. 
These effects are reversed to those of the Kerr nonlinearities which improve the 
squeezing [a]. For (5x^) « 4.44, /3 is evaluated to be about 0.88 (see Fig. [2]) 
and is in the domain given in Table [T and corresponds to the high mechanical 
displacements for the nanoresonators 17[. It also appears that for /3 = 0.1, the 
position variance {Sx'^^) is under the standard quantum limit (see Fig. [5]), allow- 
ing therefore the squeezing (up to about 31%). In order to investigate the effect 
of the optical nonlinearity rj on the mechanical squeezing, we consider rj j^ 
in the expression -^ — -^ + r] which appears in Eqs. (jlSp and P^ . Figs. [3] 
shows the effect of 77 on the position variance {Sx'^-,-^) for /? = 0.1. One remarks 
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Figure 1: Plot of the momentum variance (5p^) versus normalized detuning jr — for rj = 0, 
using experimental parameters of Ref. 15]. The value of „ — = 1 corresponds to the mo- 
mentum variance (Sp^) = 0.6362 which means that the momentum is squeezed up to about 
37%. 
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Figure 2: Plot of the position variance 5x^) versus normalized detuning ^j — for different 
values of /3 with r} = 0. The dot dashed line is plotted for /3 = 0.39 and corresponds to 
the Standard Quantum Limit (j^ = 1; {&^) = 1). The full line is plotted with f} = 0.88 

and experimental parameters of Ref. llSH and shows that the position is unsqueezed (ry- = 1; 



{Sx^) fa 4.44). The cercled line is plotted for /3 : 
{^ = 1; {Sxl,) ^ 0.69). 



: 0.1 and shows that the position is squeezed 



in Fig. [31a that the position variance {Sx^) becomes unsqueezed for high values 
of 77 [r] > 0.042). This effect of 77 is similar to that of the /3 shown in Fig. [51 It 
appears in Fig. [21b that ry shifts the optimal squeezing towards the left. Since 
77 is always present in optomechanical systems, it is then important to quantify 
it in experiments in order to determine exactly at which detuning the optimal 
squeezing can be evaluated. On the Fig. [31c where the two mentioned effects of 
rj are represented, {^x^) increases with 77 and the optimal position squeezing is 
not always at jp- = 1 but depends on the value of -q in the range yp- € [0.9; 1]. 
One also notes that the effects of 77 on the momentum variance are the same as 
these described on Figs. [31 

On the other hand, the squeezing is improved despite such nonlinearities 
in the nanobeams when the temperature is tuned down and/or the system is 
carried in a regime of strong optomechanical coupling. 



4. Squeezing of the optical output quadratures 



From the matrix equation @ we also obtain the solutions for the intracavity 
phase and intensity quadratures operators to be 



51 = -as 



and 



with 



Sip = as {GSxm + %/K(5(p'") (-in + ^) + A\/k(5/' 



03 = 



(-.o + f)= 



-in+'-] +A2 



(21) 



(22) 



(23) 



In order to analyze their squeezing, we use the well-known input-output relation 



„ out „ in I rZ ^ 
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Figure 3: Effect of optical nonlinearity 7] on the position variance for /3 = 0.1. a. Position variance versus rj shows that 
(<5x^) becomes unsqueezed when r; increases, b. shows that optimum squeezing shifts towards the left when rj increases, c. 
shows the combined effects of rj on (Sx^). 
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By using the spectral density Sa-^* {^) = ^ j^^ dwe-'(*^+")*((SA°"*(f7)(SA°"*(w)), 
we obtain the output spectrum of the intensity and phase as, 



Sio.t{a) = ai^'G'K\xeffVS^{n) + ^\n' + — 



A^ 



- aiA^K^ + A{n) - B{n) - c{n), 

^2 



S^a^t{VL) =aiG^K\ — + ^ 
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2 .,2 , /-<2 



fll 



|Xe//|^5.(f^) + y (fi^ 



— -A2 

4 



aiA^K^ + C^K^miDi^) + E{n)), 



where, 



A{n) = aAG^n„,K{nlj:j: - n^){A - n), 
B{n) = bAG'^n„iKnTeffiA - n), 



Gin) ^AG^—fl„,[2b{flljf 



D{n) = 
E{n) = 
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(2 ~^^) '*"!' 
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nc+ —d] D,ai, 
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(34) 
AK^rj , (35) 
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In Eqs. (I?71) and (I^S)) . the first terms are proportional to the position spectrum 
Sx (Eq. p4p ') and to the effective mechanical susceptibility (Eq. ((§])). These 
terms derive from the mechanical fluctuations of the oscillator. The other terms 
in (P7)) and (P5|) originate to the fluctuations of the input beam. Assuming 
that the system is in the quasi resonant regime (f2 w ^m), all the contributions 
related to the input beam fluctuations take constant values. The expressions 
(|27p and (|28p can now be integrated. By using the residues theorem and the 
Cauchy-Goursat theorem, one readily obtains 
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and 
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-{Sxl 



(39) 



where {Sx^) is given by Eq. p5|) . 

At the mechanical resonance (A « $7^) where it is estabhshed above that 
the position variance is unsqueezed ((fe^) = 4.44), we deduce from Eqs. ([55|) 
and (|39|) that both optical variances are unsqueezed. However, the effective 
detuning A w A + gnXm leads at the optical resonance (A = 0) to 
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Figure 4: Plot of the mean square fluctuation of the output light A7°"* versus the optical 
nonlinearity rj for /3 = 5.72 X lO""* . With the low input power used for the system Pi„ < 30 fiW 
[16|| . the output intensity is squeezed {|?7| < 7.44 X 10~^). 



Fig. [H shows that A/°"* is squeezed when |?7| < 7.44 x 10"^ which is ob- 
tained with the input power Pj„ < 30 (mW [16| . This means that at the optical 
resonance and for the low input power, the optomechanical cavity behaves as 
an optical filter or noise eater [lO|. Thus, the quantum fluctuations (shot noise) 
in the input coherent laser beam (A/™ = 1) are reduced after been reflected 
out of the optomechanical cavity (A/°"* < !)■ To improve this squeezing at the 
optical resonance one can increase rj (for the negative values of the detuning) 
or decrease rj (for the positive values of the detuning) (see Fig. [3|). By using 
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a°"* = ^ii^^^l (41) 



Eqs. iPc)) and ((Ml) , one can expressed output field as 



with 



A \ // A 



/K 






where e™ = \/ 'nsT'" • ^^"^ remarks that a°"* depends only on the optical 
nonlinearity which allows us to quantify his effect on the output field. Fig. [Sja 
shows that the output field decrease when 77 increases. This means that it is 
important to control the value of 77 in order to obtain the output field needed. 
This control of 77 also gives the value of the detuning at which the output field 
is optimal (see Fig. [SJb). Fig. [Sjc shows the two mentioned effects of 77 on the 
output field a°"*. 

It is also found that, as for the position, the decrease of temperature induces 
an improvement of intensity squeezing. Regarding the geometrical nonlinearity, 
it contributes to reduce the squeezing when it becomes large. But, it is generally 
weak around the optical resonance (see Table [T]) , so its effects are neglected at 
this sideband. 

Let us recall that the reduction of the nonlinearities in the optomechanical 
systems, constitutes the key of many future quantum optomechanical applica- 
tions. This requires the reduction of the mean displacement x of the nanores- 
onator or the increase of the fundamental mechanical frequency ri,„. This con- 
sists to use a very high finesse cavity which can be excited by the low input 
power [16|. 

It should be noted that squeezing of nonlinear optomechanical systems in 
which an optical cavity mode is coupled quadratically rather than linearly to 
the position of mechanical oscillator have been studied in Ref. [29| . While Sete 



and Eleuch, in Ref. 22], investigating nonlinear effects in an optomechanical 
system containing a quantum well, have found that as a result of the nonlinearity 
induced by the optomechanical coupling, the transmitted field exhibits strong 
squeezing at certain hybrid resonance frequencies and system parameters. 



5. Conclusion 

We have presented an analytical study of the geometrical and optical nonlin- 
ear effects on the optomechanical squeezing. Contrary to the Kerr nonlinearities, 
it is found that these two nonlinearities reduced the squeezing. At the detuning 
A K, f2,„ where the displacement is important, the momentum is squeezed while 
the position squeezing is very restricted for the small values of (3. This effect is 
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justified by the geometrical nonfinearity which depends on the bending moment 
of the resonator and takes into account its internal vibrations. At the optical 
resonance, the output intensity is squeezed when jryj is small (|?7| < 7.44 x 10~^). 
In a future work, we will investigate the squeezing as a function of the 
scanning frequency {Qm) in order to study the squeezing at different resonances. 
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Figure 5: Effect of optical nonlinearity rj on the output field, a. Output field versus rj shows that a°"' decreases when rj 
increases, b. The optimal output field shifts towards the left when t) increases, c. Combined effects of r] on a°"*. 



